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Abstract 

We describe a new framework of a sublinear expectation space and the 
related notions and results of distributions, independence. A new notion of 
G-distributions is introduced which generalizes our G-normal-distribution 
in the sense that mean-uncertainty can be also described. W present our 
new result of central limit theorem under sublinear expectation. This 
theorem can be also regarded as a generalization of the law of large number 
in the case of mean-uncertainty. 

1 Introduction 

The law of large numbers (LLN) and central limit theorem (CLT) are long and 
widely been known as two fundamental results in the theory of probability and 
statistics. A striking consequence of CLT is that accumulated independent and 
identically distributed random variables tends to a normal distributed random 
variable whatever is the original distribution. It is a very useful tool in finance 
since many typical financial positions are accumulations of a large number of 
small and independent risk positions. But CLT only holds in cases of model 
certainty. In this paper we are interested in CLT with mean and variance- 
uncertainty. 

Recently problems of model uncertainties in statistics, measures of risk and 
superhedging in finance motivated us to introduce, in [13] and [H] (see also 
[TTj . [T2] and references herein), a new notion of sublinear expectation, called 
"G-expectation" , and the related "G-normal distribution" (see Def. 14. 5[) from 
which we were able to define G-Brownian motion as well as the corresponding 
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stochastic calculus. The notion of G-normal distribution plays the same impor- 
tant rule in the theory of sublinear expectation as that of normal distribution 
in the classic probability theory. It is then natural and interesting to ask if 
we have the corresponding LLN and CLT under a sublinear expectation and, 
in particular, if the corresponding limit distribution of the CLT is a G-normal 
distribution. This paper gives an affirmative answer. We will prove that the 
accumulated risk positions converge 'in law' to the corresponding G-normal dis- 
tribution, which is a distribution under sublinear expectation. In a special case 
where the mean and variance uncertainty becomes zero, the G-normal distribu- 
tion becomes the classical normal distribution. Technically we introduce a new 
method to prove a CLT under a sublinear expectation space. This proof of our 
CLT is short since we borrow a deep interior estimate of fully nonlinear PDE in 
[5]. The assumptions of our CLT can be still improved. 

This paper is organized as follows: in Section 2 we describe the framework 
of a sublinear expectation space. The basic notions and results of distributions, 
independence and the related product space of sublinear will be presented in 
Section 3. In Section 4 we introduce a new notion of G-distributions which 
generalizes our G-normal-distribution in the sense that mean-uncertainty can 
be also described. Finally, in Section 5, we present our main result of CLT under 
sublinear expectation. For reader's convenience we present some basic results 
of viscosity solutions in the Appendix. 

This paper is a new and generalized version of [15] in which only variance 
uncertainty was considered for random variables instead random random vec- 
tors. Our new CLT theorem can be applied to the case where both mean- 
uncertainty and variance-uncertainty cannot be negligible. This theorem can 
be also regarded as a new generalization of LLN. We refer to [9] and [10] for the 
developments of LLN with non-additive probability measures. 

2 Basic settings 

For a given positive integer n we will denote by (x,y) the scalar product of x, 
y G R™ and by \x\ = {x.x) 1 / 2 the Euclidean norm of x. We denote by the 
collection ofnxn symmetric matrices and by S+ (n) the non negative elements 
in S(n). We observe that S(n) is an Euclidean space with the scalar product 
(P,Q) =tr[PQ]. 

In this paper we will consider the following type of spaces of sublinear ex- 
pectations: Let Q, be a given set and let H be a linear space of real functions 
defined on such that if X\, ■ ■ ■ , X n € TL then <p(Xi, • • • , X n ) € TL for each 
f € G/.Lip(K n ) where Ci.Li P (7SL n ) denotes the linear space of (local Lipschitz) 
functions <p satisfying 

Hx) - <p(y)\ <G(i + M m + \y\ m )\x-y\, Vx,yeR n , 

for some G > 0, m £ N depending on ip. 

TL is considered as a space of "random variables" . In this case we denote X — 
(Xi, • • • , X n ) € TL n . 
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Remark 2.1 In particular, if X,Y € H, then \X\, X m 6 7i are m 7i. More 
generally tp(X)^(Y) is still inTL if if, ip G Ci.Lip(K)- 

Here we use Cj.Ltp(Bi n ) in our framework only for some convenience of tech- 
nicalities. In fact our essential requirement is that TL contains all constants and, 
moreover, X £ H implies \X\ G TL. In general Ci.Li P (R") can be replaced by 
the following spaces of functions defined on R™. 

• L°°(R™): the space bounded Borel-measurable functions; 

• C(,(R n ): the space of bounded and continuous functions; 

• (R™): the space of bounded and fc-time continuously differentiable func- 
tions with bounded derivatives of all orders less than or equal to k; 

• C un if(W n ): the space of bounded and uniformly continuous functions; 

• C&.z,ip(R n ): the space of bounded and Lipschitz continuous functions. 

Definition 2.2 A sublinear expectation E on Tt is a functional E : TC i— > R 
satisfying the following properties: for all X, Y € TL, we have 

(a) Monotonicity: If X >Y then E[X] > E[Y]. 

(b) Constant preserving: E[c] = c. 

(c) Sub-additivity: E[X] - E[Y] < E[X - Y], 

(d) Positive homogeneity: E[XX] = \E[X], VA > 0. 

(In many situation (c) is also called property of self-domination). The triple 
(Q,,7i,E) is called a sublinear expectation space (compare with a probabil- 
ity space If only (c) and (d) are satisfied E is called a sublinear 
functional. 

Remark 2.3 Just as in the framework of a probability space, a sublinear expec- 
tation space can be a completed Banach space under its natural norm ||-|| =E[|-|] 
(see Jll f -Ub y ) and by using its natural capacity c(-) induced via E[\-\] (see J4H 
and But the results obtained in this paper do not need the assumption of 
the space- completion. 

Lemma 2.4 Let E be a sublinear functional defined on (0,7i), i.e., (c) and 
(d) hold for E. Then there exists a family Q of linear functional on (f2,7i) such 
that 

E[X] := sup E[X], V£ G Q. 

EGQ 

and such that, for each X G TL, there exists a E G Q such that E[X] := -E[-X"]. 
If we assume moreover that (a) holds (resp. (a), (b) hold) for E, then (a) 
also holds (resp. (a), (b) hold) for each E G Q. 
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Proof. Let Q be the family of all linear functional dominated by E, i.e., E[X] < 
E[X], for all X £ H, E S Q. We first prove that Q is non empty. For a given 
X £ H, we denote L = {aX : a £ M} which is a subspace of TC. We define 
/ : L — > M by /[aJ] = aELY], Va S R, then /[•] forms a linear functional on 
L and / <E on L. Since E[-] is sub-additive and positively homogeneous, by 
Hahn-Banach theorem (see e.g. |T9]ppl02) there exists a linear functional E 
on TC such that E — I on L and E < E on 7i. Thus -E is a linear functional 
dominated by E such that ELY] := E[X\. We now define 

E Q [X] = sup E[X]. 

EEQ 

It is clear that Eq=E. 

If (a) holds for E, then for each non negative element X € H, for each 
E e Q, ^[X] = > -E[-X] > 0, thus (a) also holds for E. If 

moreover (b) holds for E, then for each c £ R, — £J[c] = c] < E[— c] = — c 
and i?[c] < E[c] = c, we get i?[c] = c. The proof is complete. ■ 

Example 2.5 _Fbr some ip £ C;.ijj,(M), £ E H, let <p(£) be a gain value favor- 
able to a banker of a game. The banker can choose among a set of distribution 
{F(9, ^4-)},4£B(R),6>ee of a random variable £. In this situation the robust expec- 
tation of the risk for a gamblers opposite to the banker is: 

E[<p(£)] :=sup / (p(x)F(9,dx). 

3 Distributions, independence and product spaces 

We now consider the notion of the distributions of random variables under 
sublinear expectations. Let X = (X\, • • • , X n ) be a given n-dimensional random 
vector on a sublinear expectation space (f2i,Wi,E). We define a functional on 
Ci. Lip (R n ) by 

f x [<p] -=EMX)} :^a. iip (l") ^ (-00,00). (1) 

The triple (R ra , Ci.Li P (K"), Fx[']) forms a sublinear expectation space. Ex is 
called the distribution of X. 

Definition 3.1 Let X\ and X2 be two n-dimensional random vectors defined 
respectively in sublinear expectation spaces (f2i,7ii,Ei) and (f22,7^2,E2). They 

are called identically distributed, denoted by X\ = X2, if 

E 1 [ ( p(X 1 )] = E 2 [<p(X 2 )}, Vcp £ C LLip (R n ). 

It is clear that X\ = X 2 if and only if their distributions coincide. 
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Remark 3.2 If the distribution Ex of X € H is not a linear expectation, then 
X is said to have distributional uncertainty. The distribution of X has the 
following four typical parameters: 

fi:=E[X}, £-.= -R[-X], a 2 :=E[X 2 }, a 2 := -E[-X 2 }. 

The subsets [//,, p] and [cr 2 , a 2 ] characterize the mean-uncertainty and the variance- 
uncertainty of X. The problem of zero-mean uncertainty have been studied in 
[PS], [P4]- In this paper the mean uncertainty will be in our consideration. 

The following simple property is very useful in our sublinear analysis. 

Proposition 3.3 Let X,Y eH be such that E[Y] = -E[-Y], i.e. Y has no 
mean uncertainty. Then we have 

E[X + Y] = E[X]+E[Y]. 

In particular, ifE[Y] = E[-Y] = 0, then E[X + Y]= E[X]. 

Proof. It is simply because we have E[X + Y] < E[X] + E[Y] and 

E[X + Y}> E[X] - E[-Y] = E[X] + E[Y]. 

■ 

The following notion of independence plays a key role: 

Definition 3.4 In a sublinear expectation space (O, H,E) a random vector Y = 
(Yi, • • • ,Y n ), Yi G 7i is said to be independent to another random vector X = 
(Xi, • • • , X m ), X l eH under E[-] if for each test function tp G C LLl p(M. m x E") 
we have 

E[ip(X,Y)} =E[E[<p(x,Y)] x=x }. 

Remark 3.5 In the case of linear expectation, this notion of independence is 
just the classical one. It is important to note that under sublinear expectations 
the condition 'Y is independent to X " does not implies automatically that 'X 
is independent to Y". 

Example 3.6 We consider a case where X, Y G 7i are identically distributed 
and E[X] = E[-X] = but a 2 = E[X 2 } > q 2 = -E[-X 2 }. We also assume 
that E[\X\] = E[X+ + X~] > 0, thus E[X+] = ±E[\X\ + X] =±E[\X\] > 0. In 
the case where Y is independent to X , we have 

E[XY 2 } = E[X+a 2 - X-a 2 ] = (a 2 - a 2 )E[X+] > 0. 

But if X is independent to Y we have 

E[XY 2 } = 0. 
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The independence property of two random vectors X, Y involves only the 
joint distribution of (X, Y). The following construction tells us how to construct 
random vectors with given sublinear distributions and with joint independence. 

Definition 3.7 Let (Qj,Wj, Ej), i = 1,2 be two sublinear expectation spaces. 
We denote by 

Hi x H 2 := {Z(u u u2) = <p(X(u)i),Y{u2)) : (wi,w 2 ) eiiix Sh, 

(X,Y) G (Hi)" 1 x (H 2 ) n , V e C LLip (R m x R"), m,n= 1,2,- ••}, 

and, /or each random variable of the above form Z(wi,w 2 ) — <£>(X(wi), y(w 2 )), 

(Ei x E 2 )[Z] := Ei[(^(X)], w/iere p(ar) := E 2 [<p(:r, y)], a; G K m . 

It is easy to check that the triple (f2i x f2 2 ,7ii x 7i 2 , .-i x E 2 ) forms a sunlin- 
ear expectation space. We call it the product space of sublinear expectation of 
(Oi,Wi, Ei) and (0 2 ,W 2 , E 2 ). In i/iis wa?/ we can define the product space of 
sublinear expectation 

n n n 

i— 1 i— 1 z— 1 

o/ ant/ given sublinear expectation spaces (f2j,7ij, Ej), i = 1,2, • • • , n. In par- 
ticular, when (Qj, Wj, Ej) = (fli, Wi, Ei) we ftaae i/ie product space of the form 

(nf n ,wf n , Ef™). 

The following property is easy to check. 

Proposition 3.8 Lef X; oe ni- dimensional random vectors in sublinear expec- 
tation spaces (Qj, Hi, Ej), for i = 1, • • • , n, respectively. We denote 

Yi(wi, ■ ■■ ,ui n ) := -X"j(wj), i = 1, - • • ,n. 

Then Yi, i = 1, • • • , n are random variables in the product space of sublinear 

n n n 

expectation (J^Jfij, J^Wj, J^JlEj). Moreover we have Yi = Xi and Y i+1 is inde- 

i=i i=i i=i 
pendent to (Yi, ■ ■ ■ , Yi), for each i. 

Moreover, if (fij, Wi,Ej) = (Qi,Wi,Ei) and = Xl, /or a/Z i, iftera we a/so 
have y = y . In t/iis case y ia ca^ed independent copies of y /or z = 2, • • • , n. 

The situation U Y is independent to X" often appears when Y occurs after 
X, thus a very robust sublinear expectation should take the information of X 
into account. We consider the following example: Let Y = 9), ip € Cb(M. 2 ), 
where £ and X are two bounded random variables in a classical probability space 
(fi, T, P) and # is a completely unknown parameter valued in a given interval 
[a,b]. We assume that £ is independent of X under P in the classical sense. 
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On the space 7Y) with TL := {(p(X,Y) : ip £ C;.Lip(R 2 )}, we can define the 
following three robust sublincar expectations: 

E 1 [<p(X,Y)}= sup Ep[<p{X,1>{ZM M<p(X,Y)]=E P [mp ip(X,^,6)], 

8e[a,b] Oe[a,b] 

E 3 [ip(X,Y)} = E P [{ sup E P [<p(x,il>{Z,9)]}*=xl 
ee[a.b] 

But it is seen that only under the sublinear expectation E3 that Y is independent 
to X. 

Remark 3.9 It is possible that the above parameter 8 is in fact a function of X 
and £: 9 = Q(X, £) where & is a completely unknown function valued in [a,b], 
thus Y — tp(£, Q(X,£)) is dependent to X in the classical sense. But since O is 
a completely unknown function a robust expectation is E3. 

Definition 3.10 A sequence of d- dimensional random vectors {r]i} c *L 1 in Ti. is 
said to converge in distribution under E if for each if £ C(,(M n ) the sequence 

|E[<p(77i)]j converges. 

4 G-distributed random variables 

Given a pair of d-dimensional random vectors {X, Y) in a sublinear expectation 
space (Q,H,E), we can define a function 

G(p,A):=E[±(AX,X) + (p,Y)}, (p, A) e S(d) x R d (2) 

It is easy to check that G : R d x S(d) n R is a sublincar function monotonic in 
A £ S(d) in the following sense: For each p,p £ M. d and A, A £ S(d) 

G{p + p,A + A) <G{p,A) + G(p,A), 

G(Xp,XA) =XG(p,_A), VA>0, (3) 
G(p,A) >G( P ,A), if A> A. 

G is also a continuous function. 

The following property is classic. One can also check it by using Lemma 12.41 

Proposition 4.1 Let G : M. d x S(d) 1— > R be a sublinear function monotonic 
in A £ S(d) in the sense of (0j and continuous in (0,0). Then there exists a 
bounded subset 6 £ R d x R dxd such that 

G(p t A)= sup [ltr[AQQ T ] + {p, q )}, V(p, A) £ R d x S(d). 
(i,Q)e& 1 

The classical normal distribution can be characterized through the notion 
of stable distributions introduced by P. Levy [6] and [7]. The distribution of 
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a d-dimensional random vector X in a sublinear expectation space (f2,?i,E) is 
called stable if for each a, b£ R d , there exists c GR d and d € R such that 

<a,X) + (b,X)^(c,X)+a\ 

where X is an independent copy of X. 

The following G-normal distribution plays the same role as normal distribu- 
tions in the classical probability theory: 

Proposition 4.2 Let G : R d x S(d) i— > M 6e a given sublinear function mono- 
tonic in A £ §(d) ifte sense of (0) flftd continuous in (0,0). TTien i/iere exists 
a pair of d-dimensional random vectors (X, Y) in some sublinear expectation 
space (0,W,E) satisfying (0j and the following condition: 

(aX + bX, a 2 Y + b 2 Y) = (\/a 2 + IPX, (a 2 + b 2 )Y ), Va,6>0, (4) 

where (X,Y) is an independent copy of [X, Y). The distribution of {X, Y) is 
uniguely determine by G. 

Example 4.3 For the sublinear function G : R d i— > M defined by G(p) :— 
G(p,0), p £ R d , we can concretely construct a d-dimensional random vector 
Y in some sublinear expectation space (f2,7i,E) satisfying 

G(p):=E[{p,Y)], P eR d (5) 
and the following condition: 

a 2 Y + b 2 Y = (a 2 + b 2 )Y, Va, b e R, (6) 

where Y is an independent copy of Y . In fact we can take fl = R d , H = 
Gi.Lip(R d ) and Y(oj) = u>. To define the corresponding sublinear expectation E, 
we apply Proposition ^. 1\ to find a subset Q £ M. d such that 

G(p) = sup (p,q), P eR d . (7) 
gee 

Then for each !; £H of the form £(uj) = (p(oj), (p G Ci.Li P {R d )- w € R d we set 

E[£] = sup (8) 

It is easy to check that the distribution of Y satisfies 0) and (OJ). It is the 
so-called worst case distribution with respect to the subset of mean uncertainty 
O. We denote this distribution byU(Q). 

Example 4.4 For the sublinear and monotone function G : §(d) <— > R defined 
by G(A) := G(0, A), A € §(d) the d-dimensional random vector X in Proposi- 
tion \4-£\ satisfies 

6{A):=~E[{AX,X)], P eR d (9) 
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and the following condition: 

aX + bX = yja? + b 2 X, Va, b £ M, 



(10) 



where X is an independent copy of X . In particular, for each components Xj 
of X and X t of X, we have \f?E[Xi\ = EX* + X] = 2EX*] and \fm[-Xi] 
Y — Xi] = 2E[— Xj] if follows that X has no mean uncertainty: 

E[Xi] = E[-Xi] = 0, i = l, •••,(/. 

On the other hand, by Proposition ^. 1\ we can find a bounded subset £ §+(cf) 
such that 

-E[(AX, X)] = G(A) = - sup tr[AQ], A £ S(d). (11) 
2 2 Qee 

J/0 is a singleton = {Q}, t/ien X is a classical zero-mean normal distributed 
with covariance Q. In general characterizes the covariance uncertainty of X . 

Definition 4.5 (G- distribution) The pair of d-dimensional random vectors 
(X, Y) in the above proposition is called G- distributed. X is said to be G-normal 

distributed. We denote the distribution of X by X =^(0,0). 

Proposition 14.81 and Corollary 14.91 show that a G-distribution is a uniquely 
dchncd sublincar distribution on (R 2d , C;.L !p (]R 2d )). We will show that a pair 
of G-distributed random vectors is characterized, or generated, by the following 
parabolic PDE defined on [0,oo) x R d x R d : 

d t u — G(D y u, D 2 x u) = 0, (12) 

with Cauchy condition u\ t= o = if, where D y = {d Vi )f =l , J? 2 = {d^x )i,j=x- 
(TT2"|) is called the G-heat equation. 

Remark 4.6 We will use the notion of viscosity solutions to the generating 
heat equation Uty) . This notion was introduced by Crandall and Lions. For the 
existence and uniqueness of solutions and related very rich references we refer to 
Crandall, Ishii and Lions fSjj (see Appendix for the uniqueness). We note that, 
in the situation where a 2 > 0, the viscosity solution \12) becomes a classical 
C 1+ ~' 2+a -solution (see JBjl and the recent works of and /Jffl). Readers can 
understand t!2\) in the classical meaning. 

Definition 4.7 A real-valued continuous function u £ G([0,T] x M. d ) is called 
a viscosity subsolution (respectively, supersolution) of MSi) if, for each function 
^£ G b 3 ((0,oo)xM d xR d ) and for each minimum (respectively, maximum) point 
(t, x, y) £ (0, oo) x M d x M d of ' ip — u, we have 

d t ip - G(D y ip,D 2 x ip) < {respectively, > 0). 

u is called a viscosity solution of \l c 2\) if it is both super and subsolution. 
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Proposition 4.8 Let (X, Y) be G -distributed. For each <p £ C LLi p(M. d x R d ) 
we define a function 

u(t, x, y) := E[ip((x + VtX, y + tY)], (t, x) G [0, oo) x M. 

Then we have 

u(t + s,x,y)=E[u(t,x + VsX,y + sY)], s > 0. (13) 

We also have the estimates: For each T > there exist constants C, k > such 
that, for all t, s £ [0, T] and x, y £ K, 

\u(t,x,y) - u(t,x,y)\ < C(l + \x\ k + \y\ k + \x\ k + \y\ k )\x - y\ (14) 

and 

\u(t, x, y) ~ u (t + s,x + y)\ < £7(1 + |^| fe + \y\ k )(s + Is] 1 / 2 ). (15) 

Moreover, u is the unique viscosity solution, continuous in the sense of Jj^[ ) and 
{Hp, of the generating PDE fig)) . 

Proof. Since 

u(t, x, y) - u(t, x, y) = %p{x + VtX, y + tY)] - E[cp(x + VtX, y + tY)] 

< E[<p(x + VtX, y + tY)- ip{x + VtX, y + tY)] 

< E[C\(1 + \X\ k + \Y\ k + \x\ k + \y\ k + \x\ k + \y\ k )} 

x (\x - x\ + \y - y\) 

< £7(1 + \x\ k + \y\ k + \x\ k + \y\ k )(\x -x\ + \y- y\). 

We then have flip. Let {X, Y) be an independent copy of {X, Y). Since {X, Y) 
is £7-distributed, then 

u(t + s, x, y) = E[ip(x + Vt + sX, y+{t + s)Y)] 

= E[ip(x + V^X + VtX, y + sY + tY)] 

= E[E[(p{x + V~*x + VtX, y + sy + t9)] iStS)={x ,Y)] 

= E[u{t,x + V~sX,y + sY)]. 

We thus obtain (|T3j) . From this and (O it follows that 

u{t + s, x, y) — u(t, x, y) = E[u(t, x + V^X, y + sY) — u(t, x)] 
< E[Ci(l + \x\ k + \y\ k + \X\ k + |r| fc )(^l^l + s\Y\)]. 

Thus we obtain $T5§. Now, for a fixed (t,x,y) £ (0, oo) x R d x R d , let V G 
£7 b 1,3 ([0,oo) x R d x R d ) be such that -0 > it and ip(t,x,y) = u(t,x,y). By 
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and Taylor's expansion it follows that, for S £ (0,i), 

< E[ip(t -S,x + \/8X, y + SY)- ip(t, x, y)] 
<C(5 3 / 2 + S 2 )-d t ^(t,x,y)5 

+ ±[(D x iP(t, x,y),X)x/S + (D y ip(t, x,y),Y)5+^ (D 2 x ^(t, x, y)X, X) 6} 

= -d t ip(t, x, y)5 + E[(D y iP(t, x, y), Y) + 1 (D x ip(t, x, y)X, X)]8 + C(S'^ 2 + S 2 ) 

= -d t ^(t, x, y)S + SG(D y ^, D*V0 [t, x, y) + C(S^ 2 + S 2 ). 

From which it is easy to check that 

[d t iP-G(D y iP,D 2 x iP)}(t,x,y)<0. 

Thus u is a viscosity supersolution of (|12|) . Similarly we can prove that u is a 
viscosity subsolution of (fl2|) . ■ 

Corollary 4.9 // both (X, Y) and (X,Y) are G- distributed with the same G, 
i.e., 

G(p, A) := E[i (AX, X) + (p, Y)} = E[l (AX, X) + (p, Y)}, V(p, A) e S(d) xt d . 

i/ien = (X,?). In particular, X = -X. 

Proof. For each tp G C i . Li p(K' i x R d ) we set 

tt(i, ar, y) := E[^(a; + VtX, y + ty)] , 

x, y) := E[tp(x + y/tX, y + tY)}, (t, x) € [0, oo) x M. 

By the above Proposition, both u and it are viscosity solutions of the G-heat 
equation (|12p with Cauchy condition u\t—o = u\ t=0 = tp. It follows from the 
uniqueness of the viscosity solution that u = u. In particular 

E[tp(X,Y)}=E[tp(X,Y)]. 

Thus (X,Y) = (X,Y). m 

Corollary 4.10 Let(X,Y) be G -distributed. For each ip € C;.Li P (M ) we define 
a function 

v(t, x) := E[ip((x + VtX + tY)}, (t, x) £ [0, oo) x R d . 

Then v is the unique viscosity solution of the following parabolic PDE 

d t v-G{D x v,D 2 x v) = {), v\ t=0 = i>. (16) 

Moreover we have v(t,x + y) = u(t,x,y), where u is the solution of the PDE 
ilty) with initial condition u(t,x, y)\t=o = ^(x + y). 
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4.1 Proof of Proposition 14.21 

We now proceed to prove Proposition ^. 21 Let u — u v be the unique viscosity so- 
lution of the G-heat equation dHJ) with tt v | t=0 = ip, Then we take fi = R 2d , H = 
Ci.Li P (M. 2d ), uj = (x,y) £ M. 2d . The corresponding sublinear expectation E[-] is 
defined by, for each £ S TL of the form £(w) = [<p(%,y))(x,y)es? d e C[.L ip (R 2d ), 
E[£] = 7^(1,0). The monotonicity and sub-linearity of u v with respect to ip are 
known in the theory of viscosity solution. For reader's convenience we provide 
a new and simple proof in the Appendix (see Corollarv l6.4l and Corollary 16. 5p . 
The positive homogeneity of E[-] is easy to be checked. 

We now consider a pairs of d-dimensional random vectors (X, Y)(u>) = (x, y). 
We have 

E[ip.(X, Y)\ = 0), e C LLip (R 2d ). 

In particular, just set <fio(x, y) = \ (Ax, x) + (p, y), we can check that 

{t, x,y):= G(p, A)t + - {Ax, x) + (p, y) . 

We thus have 

E[i (AX, X) + (p, Y)] = u*°(t, 0)|*=i = G(p, A), (p, A) e M d x S(n). 

To prove that the distribution of (X,Y) satisfies condition (13} , we follow 
Definition 13.71 to construct a product space of sublinear expectation 

(Q, H,E) = (fixO,HxM,ExE) 

and introduce two pair of random vectors 

(X,Y)(fJi,w 2 ) = wi, (X,Y)(w 1 ,u 2 )=U2, Kw 2 )eflxfl. 

By Proposition [3U both (X,Y) = (X,Y) and (X,Y) is an independent copy of 
{X,Y). For each ip e C LLlp (R 2d ) and for each fixed A > 0, (x,y) € R 2d , since 
the function v defined by v(t, x, y) := u v (Xt, x + y/Xx, y + Ay) solves exactly the 
same equation (fT2"j) but with Cauchy condition 

v\ t =o = <p(x + Vx x •, y + A x •). 

Thus 

E[<p(x + VXX, y + XY)} = v(t, x, y)|t=i 

= u^*-^\t,x,y)\ t=1 =v?{X,x,y). 

By the definition of E, for each t > and s > 0, 

E[tp(y/tX + yfsX, tY + sY)} = E[E[<p(Vix + yfsX, ty + sY)} {x>y)={XiY) ] 

= u uV ( s '-'-\t, 0, 0) = u v (t + s, 0, 0) 

= t{ip(VTT7sX, (t + s)Y)}. 
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Namely (y/lX + y/aX, tY + sY) = (y/TTsX, (t + s)Y). Thus the distribution 
of (X, Y) satisfies condition . 

It remains to check that the functional E[«] : Ci.n p (R 2d ) i— > K forms a 
sublinear expectation, i.e., (a)-(d) of Definition 12.21 are satisfied. Indeed, (a) 
is simply the consequence of comparison theorem, or the maximum principle of 
viscosity solution (see [CIL] , the prove of this comparison theorem as well as the 
sub-additivity (c) are given in the Appendix of [P6]). It is also easy to check 
that, when tp = c, the unique solution of (| 12(1 is also u = c; hence (b) holds true, 
(d) also holds since u Xv = Xu v , A > 0. The proof is complete. 



5 Central Limit Theorem 

Theorem 5.1 (Central Limit Theorem) Let a sequence {(Xi,Yi)}'^_ 1 ofR d x 

R d -valued random variables in (7i, E). We assumed that (X;_|_i, li+i) = (Xj, Yi) 
and (Xi + i, Yi + i) is independent to {(Xi, Y\), • • • , (Xj, Yi)} for each i = 1,2, • • • . 
We assume furthermore that, 

E[Xi] = E[-Xi] = 0, 

Then the sequence {S n }^ =1 defined by 

^ yn n 

4 — 1 V 

converges in law to £ + £: 

lim Eb(5 n )]=E[^ + C)], (17) 

n — >oo 

/or aZ/ functions (f € G(R d ) satisfying a polynomial growth condition, where 
(£i C) * s a pa*r of G -distributed random vectors and where the sublinear function 
G : S(d) x M. d i— > R is denned 6y 

G(p,A):=E[<p,Yi) + i(Aari,Xi)], AeS(d), P el d . 

Corollary 5.2 T/ie sum X)"=i ^ converges in law toAf(0, 9), where the subset 

9 C §+(d) is defined in {77]) /or G(A) = G(_0, A), A € S(d). 77ie s«m £_™=i if 
converges in law to U(Q), where the subset cM. d is defined in /or G(p) = 
G(p, 0), p € R d . J/ we iafce in particular ip(y) — dg(y) = mf{|x — y\ : x £ ©}, 
i/ien 6?/ (0) we Ziaue the following generalized law of large number: 

n y 

lim E[de(V — )] = supdgW = °- ( 18 ) 
™^°° i=1 n see 
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Remark 5.3 If Yi has no mean-uncertainty, or in other words, is a singleton: 
= {9} then U8\) becomes 

n 

lim E[| V — - 9\] = 0. 

n — >oo ^ — ' n 

8=1 

To our knowledge, the law of large numbers with non-additive probability mea- 
sures have been investigated with a quite different framework and approach from 
ours (see fMj, fWj/). 

To prove this theorem we first give 

Lemma 5.4 We assume the same condition as Theorem 15.11 We assume 
furthermore that there exists f3 > such that, for each A, A £ S(d) with A > A, 
we have 

E[(AX 1 ,X 1 )]-E[(AX 1 ,X 1 )]>pti{A-I]. (19) 

Then (Tfy holds. 

Proof. We first prove ([TTJ) for ip £ C b . Lip (R d ). For a small but fixed h > 0, let 
V be the unique viscosity solution of 

dtV + G{DV, D 2 V) = 0, (t, x) G [0, 1 + h] x R d , V|t=i +/[ = (p. (20) 

Since (£, Q is G-distributed we have 

V(M)=Efa(£ + C)]> + = (21) 

Since ([20]l is a uniformly parabolic PDE and G is a convex function thus, by 
the interior regularity of V (see Krylov [5], Example 6.1.8 and Theorem 6.2.3), 
we have 

\\V\\ cl+ct/2 ,2 +a(l0:1]xRd) < oo, for some a £ (0, 1). 
We set S = A and 5 = 0. Then 

n-1 

S„) - V(0, 0) = + 1)8, S i+1 ) - V(iS, Si)} 

i=0 

n-1 

= J2 + w - v( - i6 > + w 5 > ^f+i) - &)]> 

i=0 
n-1 

with, by Taylor's expansion, 

J\ = d t V{i5, S t )5+^ (D 2 V{i5, Si)X i+1 ,X i+1 ) 5+^DV{i5, Si), X i+1 VS + Y +1 5^ 
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P s = I [d t V((i + (3)8, S l+1 ) - d t V(iS, S l+1 )]df36 + [dtV{iS, S i+1 ) - d t V(iS, S^S 



+ \ (D 2 V(i5, S l )X l+l ,Y l+l ) 5 Z ' 2 + l - (D 2 V(i5, S t )Y l+1 , Y i+1 ) 5 



2 

r-l pi 



^0 

with _ 

6^ 7 = D 2 V(i5, S, + 7 /3(X i+1 V5 + Y i+1 S) - D 2 V(i5, S t ). 

Thus 

71 — 1 71—1 71 — 1 71—1 

i=0 i=0 i=0 i=0 

We now prove that E^™^ 1 "^J] = ^- F° r the 3rd term of J] we have: 
E[(w(tMi),*»+iV^)] =±[-(DV(i6,Si),X i+1 Vs}] =0. 
For the second term, we have, from the definition of the function G, 
E[J|] = E[9iF(0, ft) + G(DV(iS, Si),D 2 V(iS, St))]5. 

We then combine the above two equalities with dtV + G(DV, D 2 V) — as well 
as the independence of (Xi+i, Yi+i) to {(X\, Yy), ■ • • , (Xi, !»)}, it follows that 

n— 1 n — 2 



^^ = ^^ = •■■ = 0. 

i=0 i=0 

Thus (|2"2"]) can be rewritten as 

71—1 TO — 1 

-E[- £ I}} < S„)] - 7(0, 0) < E[£ 

i=0 i=0 

But since both dtV and D 2 V are uniformly a- holder continuous in x and ~ 
holder continuous in t on [0, 1] x R, we then have 

< C5 1+a ' 2 [l + \X t+1 \ 2+a + |yi| 2+Q ]. 

It follows that 

< C5 1+a / 2 (l+E[\X 1 \ 2+a + \Y 1 \ 2+a ]). 

Thus 

-C(-)^ 2 (l + E[\X 1 \ 2+a + |Yi| 2+Q ]) < E[V(l,S n )} - V(0,Q) 
n 

< C(-) Q/2 (1+E[|X 1 | 2+Q + |Yi| 2+a ]). 
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As n 



oo we have 



lim E[V(l,S n )} =V(0,0). 

n — >oo 

On the other hand, for each t, t' £ [0, 1 + h] and x £ K d , 



we have 



(23) 



\V(t, x) - V(t', x)\ < C(y/\t=¥\. + \t-t'\). 



Thus |V(0,0) - V(h,0)\ < C(Vh + h) and, by (J23]), 

|E[y(l,5„)]-E[«p(5„)]| = \E[V(l,S n )]-E[V(l + h,S n )]\<C(Vh + h). 
It follows form ((HJ and ([23]) that 

limsup|E[p(S n )]-Efo>(£ + C)]| < 2C(Vft + /i). 



We now give 

Proof of Theorem 15.11 For the case when the uniform Elliptic condition [TO] 
does not hold, we first follow an idea of Song [17 to introduce a perturbation 
to prove the above convergence for ip £ Cb.Li P (^- d )- According to Definition 13. 71 
and Proposition 13.81 we can construct a sublinear expectation space (f2,W, E) 
and a sequence of three random vectors {(Xj, Yi, fji)}i^i such that, for each n = 



1,2, ■ {(Xi,Yi)}? =1 = {(Xi,Yi)}? =1 and (X n+1 ,Y n+1 ,r) n+1 ) is independent to 
{(X^Y^fji)}^ and, moreover, 



We then use the following perturbation Xf = Xi + er\i for a fixed e > 0. It is 
seen that the sequence {{Xf ,Y i )}°^ 1 satisfies all conditions in the above CLT, 
in particular 



Since h can be arbitrarily small we thus have 



lim E[tp(S n )] =E[y>(0]. 




G £ (p,A) := E[- (AXl,Xl) + (p,^)} = G{p,A) + -tr[A}. 



Thus it is strictly elliptic. We then can apply Lemma [5^41 to 




n 



n 



and obtain 



lim Efc>(S*)] 



E[ v (f + C + «;)]. 



n 
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where (£, () is G-distributcs and 

EMZ + C,v)] = -i= [ n^ + C,x)]e^dx, 4>eC LLvp {R 2d ). 

Thus (£ + en, () is G e -distributed. But we have 

\E[ip(S n )\ - E[p(S£)]| = |Eb(5„)] - Efo>(S„ + e J n )]| 
< eCE[\J n \] < Ce 

and similarly |E[y>(£)] — E[y>(£ + err?)] j < Ce. Sine e can be arbitrarily small, it 
follows that 

lim &[<p(S n )] = + C)}, V<P G C b . Lip (R d ). 

n — >oo 

On the other hand, it is easy to check that sup„ E[|5„|] +E[|£ + £|] < oo. We 
then can apply the following lemma to prove that the above convergence holds 
for the case where ip in C(R d ) with a polynomial growth condition. The proof 
is complete. 

Lemma 5.5 Let (£l,Tt,E) and (Cl,H,E) be two sublinear expectation space and 
letC, gH and £ H, n = 1, 2, • • • , be given. We assume that, for a given p>l 
we have sup n E[|F„|P] + E[\Y\p] < C. If the convergence lim„^oo E[ip{Y n )} = 
E[ip(Y)] holds for each if € Cb.Lip(^ d ), then it also holds for all function if G 
C(R d ) with growth condition \ip(x)\ < C(l + \x\p~ 1 ). 

Proof. We first prove that the above convergence holds for ip e C b (R d ) with 
a compact support. In this case, for each e > 0, we can find a (p <G Cb.Li P (R d ) 
such that sup^gRd \<p{x) — <p(x)\ < §• We have 

\EMY n )]-E[p(Y)}\ < \E[ip(Y n )} + \E[p(Y)}-E[p(Y)}\ 

+ \E[p(Y n )} - E[p(Y n )}\ < e + |%(y„)] - E[p(Y)}\. 

Thus limsup n _ >oc |E[(p(F„)] — E[<^(y)]| < e. The convergence must hold since e 
can be arbitrarily small. 

Now let ip be an arbitrary Cb(R")-function. For each N > we can find 
Pi,f2 € Cb(Mr) such that ip = if\ + pi where if\ has a compact support and 
ipi{x) = for \x\ < N, and |v?2(a0| < \<f( x )\ f° r au x - It is clear that 

|V2(^)| < — — — — — — - , Vx, where (7= sup 

Thus 

\E[p(Y n )]-E[p(Y)}\ = \E[ ( p 1 (Y n ) + MY n )]-n'Pi(Y) + 'P2(Y)]\ 

< |E[^(r„)] - Efo>i(y)]| + E[|^ 2 (r„)|] + e[|<^ 2 (f)|] 

< |E[^(y„)] - E[ Vl (y)]| + ^(E[|y n | +I[|y|]) 

<|E[^ 1 (y„)]-E^ 1 (y)]| + — . 
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We thus have limsup^^ |E[^(y, l )]-E[^(F)] | < Since N can be arbitrarily 
large thus E[y>(i^)] must converge to E[y(y)]. ■ 

6 Appendix: some basic results of viscosity so- 
lutions 

We will use the following well-known result in viscosity solution theory (see 
Theorem 8.3 of Crandall Ishii and Lions [I])- 

Theorem 6.1 Let Ui £ USC((0, T) x Qi) for i = 1, • • • ,k where Qi is a locally 
compact subset of M. Ni . Let cp be defined on an open neighborhood of (0, T) x 
Qi x ■ ■ ■ x Qk and such that (t, X\, ■ ■ ■ , Xk) is once continuously differentiable in 
t and twice continuously differentiable in (xi, ■ ■ ■ ,Xk) € Qi X • • • X Qk- Suppose 
that t G (0, T), Xi £ Qi for i = 1, • • • , k and 

w(t,xi,- ■■ ,x k ) := x\) -\ h u k {t,x k ) - tp(t,xi, ■■ ■ ,x k ) 

< w(i,xi, ■■ ■ ,x k ) 

for t £ (0,T) and Xi £ Qi. Assume, moreover that there is an r > such that 
for every M > there is a C such that for i = 1, • • • , k, 



bi < C whenever (bi,qi,Xi) £ V 2 ' + Ui(t, Xi), 

\xi-Xi\ + \t-i\<r and \ui(t,Xi)\ + \q, t \ + < M. 

Then for each e > 0, there are Xi £ S(Ni) such that 

2 + 

(i) (bi,D Xi ip(t,xi,- ■ ■ ,Xk), Xi) £ V Ui(t,xi), i=l, ■■■,k; 
(H) 



(24) 



(- + \\M)< 



(Hi) bi H h b k = d t if{i, xi, 

where A = D 2 tp(x) £ S(Ni + 





,Xk) 







Xk 



< A + sA 



Observe that the above conditions (|24|) will be guaranteed by having Ui be 
subsolutions of parabolic equations given in the following two theorems, which 
is an improved version of the one in the Appendix of [16] . 

Theorem 6.2 (Domination Theorem) Ui £USC([0, T] xR^) be subsolutions of 



dtu — Gi(t, x, u, Du, D 2 u) = 0, i = 1, 



(25) 



on (0, T) xl w such that, for given constants > 0, i = 1, • • • , k, (Y^—i Ui(t, x)J 

0, uniformly as \x\ — > oo. We assume that 
(i) The functions 



d : [0, T] x R N x R N 



x S(N) 



18 



are continuous in the following sense: for each t € [0,T), v\, v-i € K, x, y, p, 
qeR N andY e S(N), 

[Gi(t,x,v,p,X) - Gi(t,y,v,p,X)}~ 

< w(l + (T — ty 1 + \x\ + |y| + \v\)u(\x -y\ + \p\ ■ \x - y\) 

where u, uj : R + i— > R + are given continuous functions with lu(0) = 0. 

(ii) Given constants > 0, i = 1, • • • , k, the following domination condition 

holds for Gi'. 

k 

^>2 (3iGi(t, x,Vi,pi, Xi) < 0, (26) 

i=l 

for each (t,x) € (0, T) x R N and (v l ,p l ,X l ) e M x R N x S(iV) audi i/iaf 

Ei=i A«< > 0> Eti Aft = 0, Eti < o- 

XTien a similar domination also holds for the solutions: If the sum of initial 

values Ei=i ') * s a non-positive and continuous function on R N , then 

E^=iA«i(*r)<0>/or a/it>0. 

Proof. We first observe that for 5 > and for each 1 < i < k, the functions 
defined by Ui := Ui — S/(T — t) is a subsolution of: 

d t Ui-Gi(t,x,Ui + 8/(T ~t),Dui,D 2 Ui) < -— — 

(T-ty 



where Gi(t, x, v,p, X) := Gi(t,x,v + S/(T — t),p,X). It is easy to check that 
the functions Gi satisfy 

from Ei=2 A^i < in t 
additional assumptions 



the functions Gi satisfy the same conditions as Gi. Since J2i=i ftiUi < follows 
from Ei=2 Pi&i < in the limit 5 J, 0, it suffices to prove the theorem under the 



d t Ui - Gi(t, x, Ui, Dui, D 2 Ui) < -c, c := S/T 2 , 
and ]im t _>x Vn(t,x) = -co, uniformly in [0, T) x 

To prove the theorem, we assume to the contrary that 

ft 

sup > , Pi u i(t, x ) — m o > 
(s,i)£[0,T)xt» i=1 

We will apply Theorem lO for x = (xi, ■ ■ ■ ,x k ) eR kxN and 

k fe-1 

w(t,x) := }^iUi(t,Xi), ip{x) = tp a {x) := -^y \x t+1 - x t 
i=l i=l 



piV 



(27) 
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For each large a > the maximum of u> — ip a achieved at some (t a , x a ) inside 
a compact subset of [0, T) x R fexAr . Indeed, since 

ft 

M a =Y J PMt a ,xf) - <p a (x a ) > m , 

i=l 
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thus t a must be inside an interval [0,Tq], Tq < T and x a must be inside a 
compact set {x £ xR kxN : snp t£ ^ QTg ^w(t,x) > We can check that (see 

|2 a Lemma 3.1) 



(i) lim Q ^oo (f a (x a ) = 0. 

(ii) lim Q ^oo M a = lim^oo piui(t a ,x? ) + ■ • ■ + /3 k u k (t a , x%)) 

= su P(M)e[0,T)xR«[/3lWl(M) + ' ' ' +PkU k (t,x)] 
= [[3iui(t, x) H h f3 k u k (i, x)} = m . 



(28) 



where (i, x) is a limit point of (t a , x"). Since iti £ USC, for sufficiently large a, 
we have 

/W^r') + --- + /W^X) > 

If t = 0, we have limsup^^ J2i=i 0i u i{t a i x ?) = J2i=i ft"i(0,x) < 0. We 
know that i > and thus t a must be strictly positive for large a. It follows 
from Theorem 16 . II that . for each e > there exists bf £ M, X, £ §(7V) such that 

k 

(bf,D Xi <p(x a ), Xi) £ J 2 Q + ( Ul )(t a , xf), ^ = 0, for i = 1, • • • , fe, (29) 

i=l 

and such that 



{- + \\A\\)I< 



( &iX x 


V o 



\ 



Pk-iX k -i 

o p k x k ) 



< A + eA 



(30) 



where A = D 2 ip a (x a ) £ S(kN) is explicitly given by 



/ / 



A = aJ kN , where J kN 



N 



-In 



-In \ 



0...-I N I N -I N 
V -In ... -I N I N J 



The second inequality of (I30|) implies J2i=i PiXi < 0. Setting 
P? = A^M*") = /3r 1 a(2 a; ? - x£ - a£), 

Pfc - D Xh ip a (x a ) = (3 k 1 a(2x a k - x%_ x - <). 
Thus X)* =1 &P? = 0. This with $29]) and JS7J) it follows that 



6f-G i (t°,a:?,« i (t 0, ,a ; f),p?,X i )<-c i = !,•■■,*;. 
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By ([2"8])-(i) we also have linia^oo \pf\ ■ \xf — x"\ — > 0. This, together with the 
the domination condition (|26[) of Gi , implies 

k k 

k 

> -^2(3iGi(t a ,x",Ui(t a ,Xi),pf,Xi) 

i=i 

k 

-Y,^[G l {t a ,xt,u l {t a ,x':), P t,X. l )~G % {t a ,x'i,u l {t a ,xr),p^X l )]- 

i=l 
k 

> -£>w(l + (T - To)- 1 + |^?| + |<| + Itii^.OIJwda:? - s?| + |pf | • |< - 

s=i 

The right side tends to zero as a — > oo, which induces a contradiction. The 
proof is complete. ■ 

Theorem 6.3 (Domination Theorem) Let polynomial growth functions m £USC([0,T] 
M. N ) be subsolutions of 

d t u-Gi(u,Du,D 2 u) = 0, i=l, (31) 

on (0, T) xR N . We assume that G, : R x M. N x §(N) >->■ K, i = 1, • • • , k, are 
given continuous functions satisfying the following conditions: There exists a 
positive constant C , such that 

(i) for all A > 

XGi(v,p,Y) < CG t (Xv,Xp,XY); 

(ii) Lipschitz condition: 

\Gi( Vl ,p,X) - Gi(v 2 ,q,Y)\ < C(\vi - v 2 \ + \p - q\ + \\X - Y\\), 
Vvi,V2 6R, Vp,q£R N and X,Y £§(N), 

(Hi) domination condition for G{: for fixed constants fa > 0, i = 1, • • • , k, 

k 

^fcGiivuPuXi) < 0, for all v t £l, K 6 R N , X, e S(JV), 
»=i 

such that @i v i ^ °> fcP* = °' @ iXi - °' 

i— 1 i— 1 i— 1 

Then the following domination holds: If X)j=i A'Ui(0, •) is a non-positive 
continuous function, then we have 

k 

Y, /W*> x) < 0, V(t, x) e (0, T) x R N . 
i=i 
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Proof. Wc set £0) := (1 + \x\ 2 ) 1 / 2 and 

Ui(t, x) := Ui(t, x)e xt £,~ 1 (x), i = 1, ■ ■ ■ , k, 

where I is chosen large enough so that 53i=i l^itA^OI — > uniformly. From 
condition (i) it is easy to check that for each i = 1, • • • , Ui is a subsolution of 

d t Ui - G;(a;,i2;,L>i2;,Z) 2 u;) = 0, (32) 

where 

Gi(x,v,p,X) := -Xv 

+ e xt CG l {v,p + vn(x),X + p <g> n(x) + n(x) ® p + vk(x)). 

Here 

T]{x) := C 1 {x)D^(x) = k(l + \x\ 2 )~ l x, 

k(x) := C 1 ^)!) 2 ^) = ^(1 + M 2 )^ - fc(fc - 2)(1 + |a;| 2 ) _2 a; ® a;. 

Since 77 and k are uniformly bounded, one can choose a hxed but large enough 
A > such that Gi(x, v,p, X) satisfies all conditions of G{, i = 1, ••• ,k in 
Theorem 16.21 The proof is complete by directly applying this theorem. ■ 
We have the following Corollaries which are basic in this paper: 

Corollary 6.4 (Comparison Theorem) Let F\,F2 : K w x E>(N) i— » K 6e given 
functions satisfying similar conditions (i) and (ii) of Theorem \ 6.3[ We also 
assume that, for each p,q £ M w and X, Y £ S(iV) suc/i £/ia£ 1 > 7, me /iawe 

Fi(p,X)> F 2 (p,y). 

Lei £LSC((0,T) x l 1 *), i = 1,2, &e respectively a viscosity supersolution of 
dtv — Fi{Dv, D 2 v) = smc/i i/iai «i(0, •) — ^2(0, •) is a non-negative continuous 
function. Then we have v\{t,x) — V2(t,x) > /or a/Z (t, x) £ [0, oo) x $L N . 

Proof. We set [3 X = j3 2 = I, Gi(p,X) := -Fi(-p,-X) and G 2 = F 2 {p,X). 
It is observed that u\ := —vi eUSC((0,T) x R N ) is a viscosity subsolution of 
d t u - Gi(Du,D 2 w) = 0. For each p u p 2 £ R N and X u X 2 £ S(N) such that 
Pi + p 2 = and X\ + X 2 < 0, we also have 

G 1 ( Pl ,X 1 ) + G 2 (p 2 ,X 2 )=F 2 (p 2 ,X 2 )-F 1 (p 2 ,~X 1 )<0 

We thus can apply Theorem l6.3l to get —ui + u 2 < 0. The proof is complete. ■ 

Corollary 6.5 (Domination Theorem) Let F, : R N x S{N) i-> R, i = 0,1, 
given functions satisfying similar conditions (i) and (ii) of Theorem \ 6.S\ Let 
Vi £LSC((0,T) x R N ) be viscosity super solutions of dtv — Fi(Dv, D 2 v) = 
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respectively for i = 0, 1 and let v 2 G USC((0, T) x R N ) be a viscosity subsolution 
of d t v — Fi(Dv, D 2 v) = 0. We assume that: 

F 1 {p,X)-F 1 {q,Y)<F {p-q,Z), 

Vp, q G R N , X,Y,Z G S{N) such that X -Y < Z. 

Then the following domination holds: If vq(0, •)+i?x(0, ■) — ^2 (0, •) is a continuous 
and non-negative function then Vo(t, •) +v\(t, •) — v 2 (t,-) > for all t > 0. 

Proof. We denote 

G l {p,X):=-F l (-p,-X), * = 0,1, and G 2 (p,X) := FfaX). 

Observe that Uj = — Wj GUSC((0,T) x R*), i = 0, 1, are viscosity subsolutions 
of <9 t u - Gi(£>u, L> 2 u) = 0, i = 0, 1. We thus have, for each X + X l + X 2 < 0, 
Pa + Pi + P2 = 0, 

G (p , X ) + dipuXx) + G 2 ( P2 , X 2 ) 
= -F {-p ,-X ) -F 1 (- Pl ,-X 1 ) +F 1 (p 2 ,X 2 ) < 0. 

Theorem l6.3l can be applied, for the case Pi = 1, to get ^ Ui < 0, or vo+vi— v 2 > 
0. ■ 

Another co-product of Theorem 16.31 is: 

Corollary 6.6 (Concavity) Let F : R N x S(N) ^ R be a given function sat- 
isfying similar conditions (i) and (ii) of Theorem 1 6. 31 We assume that F is 
monotone in X, i.e. F(p,X) > F(p,Y) if X > Y , and that F is concave: for 
each fixed a G (0, 1), 

aF(p, X)+(l-a)F(q, Y) < F{ap+{l-a)q, aX+(l-a)Y), Vp, q G R N , X, Y G S(N). 

Let Vi EUSC((0,T) x R N ), i — 0,1, be respectively viscosity subsolutions of 
dtv — F(Dv, D 2 v) = and let v ELSC((0, T) x R n ) be viscosity supersolution of 
dtV — F{Dv, D 2 v) = such that avi(0, -) + (l — a)v 2 (0, ■) — v(0, •) is a non-positive 
continuous function. Then for all t > av%(t, •) + (1 — a)v 2 (t, ■) — v(t, •) > 0. 

Proof. We set 0\ = a, 02 = (1 — a), 0s = 1 and denote 

Gi(p,X) = G 2 (p,X) := F(p,X), G 3 (p,X) = -F(-p,-X). 

Observe that m = Vi GlISC((0,T) x R n ), i = 1,2, are viscosity subsolutions 
of d t u — Gi(Du,D 2 u) = 0, M3 = — v G USC is a viscosity subsolution of d t u — 
Gs(Du^ D 2 u) = 0. Since F is concave, thus for each Pi G 1* and Xi G S(JV) 
such that P1X1 + [3 2 X 2 + f3 3 X 3 < 0, fiipi + /3 2 p 2 + /3 3 p 3 = 0, we have 

f3 1 G 1 (p 1 ,X 1 ) +p 2 G 2 (p 2 ,X 2 ) +p 3 G 3 (p 3 ,X 3 ) < F(f3 lPl + fop 2 , ^Xt + p 2 X 2 ) - F{-p 3 

< F(- P3 , [3 1 X 1 + (3 2 X 2 ) - F(- P3 , -X 3 ) 

< 0. 

Theorem [O] can be applied to prove that av\(t, •) + (1 — a)v 2 (t 7 •) < v(t, •). The 
proof is complete. ■ 
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